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A SIMPLIFIED METHOD FOR CALCULATING THE ATMOSPHERIC HEATING RATE BY 


ABSORPTION OF SOLAR RADIATION IN THE 
STRATOSPHERE AND MESOSPHERE 
Tatsuo Shimazaki and Leland C. Helmle* 
Ames Research Center 


SUMMARY 


A simple analytical formula is worked out for calculations of the 
atmospheric heating rate by absorption of solar radiation by O 3 , H 2 O, and 
CO 2 . The method needs only seven parameters for each molecule and can save 
computational time and memory locations. It is particularly useful for 
heating calculations in three-dimensional global circulation models below 
80 km. Applying the formula to the observed distributions of O 3 , H 2 O, and 
CO 2 produces reasonable latitudinal and seasonal variations in the heating 
rate. The calculated heating rate, however, is sensitive to the global dis- 
tributions of the absorbing gases, and uncertainties in the O 3 distribution 
above ^50 km and the H 2 0 distribution below ~20 km may seriously affect the 
global distributions of the heating rate in these regions. 


INTRODUCTION 


Heating the atmosphere by solar radiation of wavelengths shorter than 
~4 y is the ultimate source of energy to drive the atmospheric motion. About 
99 % of the solar radiation energy is in this spectrum range with the maximum 
radiation at ~4828 A corresponding to the 6000 K black body radiation. The 
main radiative processes in the earth's atmosphere for this spectrum range are 
absorption, reflection, and scattering; emission is entirely negligible 
because the atmospheric temperature is too low. 

Calculations of the atmospheric heating rate due to gaseous absorption 
of solar radiation are essential in studies of radiative balance and global 
circulation models of the atmosphere. The main absorbing molecules are water 
vapor (H 2 O) in the troposphere and ozone (O 3 ) in the stratosphere and meso- 
sphere, whereas carbon dioxide (C0 2 ) is of secondary importance as an absorber 
throughout the region from troposphere to mesosphere. Molecular oxygen (0 2 ) 
becomes the dominant absorber in the thermosphere. The absorption of each 
constituent is highly selective and strongly dependent on wavelength, and cal- 
culations of the heating rate taking into account the detailed wavelength 
dependence are very time consuming if not impossible. To avoid such incon- 
veniences and to facilitate the heating rate calculations, empirical formulae 
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have been developed by various investigators. Formulae differ widely, 
however, and in most cases different formulae have been used for different 
constituents and for different wavelength ranges even by the same author (s). 
Thus, we need different computer programs for calculating the heating rate in 
each case. 

The purposes of this technical paper are (1) to calculate the heating 
rate due to O 3 absorption with high precision and to formulate this result 
and experimental results of heating due to H 2 0 and C0 2 absorption by 
Howard et al. (refs. 1 and 2) into a simple analytical expression; (2) to 
survey the existing formulae for heating rate calculations and compare our 
result with the results calculated by these formulae; and (3) to illustrate 
the global distribution of the heating rate calculated by our formula using 
the global data on O 3 , H 2 0, and C0 2 density distributions . Development of 
such an analytical formula is useful for simplifying the radiation calcula- 
tions in global circulation models, since the computation time and storage 
locations are very limited by the fact that a large amount must be reserved 
for calculations of dynamics. 

We intend to use the same formula used in calculations of the photodisso- 
ciation rates of molecules at the Schumann-Rung e band (Shimazaki et al., 
ref. 3), so that the heating rate and the photodissociation rate can be cal- 
culated by the identical formula. We need only several parameters in each 
calculation, and the same short subroutine can be used for calculations of 
both the heating rate and the photodissociation rate, thereby saving both 
computation time and storage locations. 


METHOD OF ANALYSES 


Heating by Ozone 

Ozone absorbs solar radiation energy mainly at the Hartley band 
(2400 to 3000 A), Huggins band (3000 to 3600 A) and Chappuis band 
(4000 to 8500 A) „ The division at 3000 A between the Hartley and Huggins 
bands is somewhat arbitrary, and the Hartley band includes the wavelengths 
shorter than 2400 A. The absorption by 0 2 becomes increasingly important 
below 2400 A, but heating by those absorptions should be small in the atmo- 
spheric region under our present concern. Since the absorption cross sec- 
tions of ozone are well determined for the entire spectrum range, the exact 
heating rate can be calculated straightforwardly by 

Q(X,u) = I(X)a(X)e _a(A)u (1) 

where A is the wavelength, u the total amount of ozone in cm atm NTP , I 
the incident solar energy flux, and a the absorption cross section. I is 
related to the photon flux <f> by the relation 

I(X) = ^ * (2) 
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where h is the Planck’s constant and c the speed of light * In the region 
of interest (below 80 km), the 0( 3 P) formed by the ozone photolysis quickly 
recombines; the portion of energy spent for dissociation will be regained for 
atmospheric heating through chemical recombination processes. Thus, the 
energy loss by dissociation can be neglected. 

The absorbed energy flux along the path is then calculated by 

S(A,u) = f U Q(A,u)du = 1(A) [1 - e " CT(X)u ] (3) 

•'o 

In calculating equation (3) we used data for the solar flux from Detwiler 
et al. (ref. 4) and Ackerman (ref. 5) for the 2400 to 3000 A range and from 
Thekaekara and Drummond (ref. 6) for the 3000 to 8500 A range, and the data 
for the absorption cross section from Inn and Tanaka (ref. 7). The entire 
spectrum range 2400 to 8500 A is divided into 122 subdivisions with the equal 
interval (50 A) , and the integrated solar flux and the averaged cross section 
have been calculated for each subdivision. 

Taking the summation of equation (3) for the pertinent spectrum range 
would give the absorbed energy at the respective band 

S(u) =£ S(A,u) (4) 

A 

The calculated absorbed energy is shown in figure 1 for the Hartley band, 
Huggins band, and Chappuis band, as well as for the total of the three bands. 
The abscissa indicates the ozone column density measured in cm atm NTP, and 
the corresponding scale for the height (km) is shown at the top for the three 
representative seasons and latitudes. In calculating these height scales, 
the diurnally averaged mass factor is multiplied by the vertical column den- 
sity to obtain the total ozone density in the column along the oblique path 
in each case. Effects of the dependence of absorption on pressure are also 
incorporated for evaluating these height scales (see the later discussion on 
the effective mass for the pressure dependence of absorption) . 

It is seen in figure 1 that the solar energy is absorbed mainly at the 
Hartley band for u less than ~0.1 cm, whereas the absorption occurs mainly 
at the Chappuis band for u greater than ~1 cm. In the transition region 
there is a depression for the total of the three bands. The result of the 
calculations shown in figure 1 is fitted to the analytical expression 

o k 

log S = C 0 + Cj log u + C£ (log u) 2 + . . .4- C (log u) (5) 

K, 

for the range of u between 10” 7 and 10 2 cm. The coefficients are 

determined by the least square method so that the variations of the curves 
are best fitted with the equation (5) . The resultant coefficients are 
tabulated in table 1(a) for the case of k = 7; for the total of three bands 
the result calculated from equation (5) is shown in figure 1 by a dash-dotted 
curve, which agrees well with the original solid curve. For the convenience 
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of some modelers, coefficients determined for the MKS system of units are 
given in table 1(b). The mean error of log S, or the mean relative error of 
S, is defined by 

E(k) = ( lo S Slf ~ log s 5 ) 2 /( n ~ k ) ( 6 ) 

where S 4 and S 5 are the values of S calculated by equations (4) and (5), 
respectively, n is the total number of observational points and k is the 
degree of the polynomial. 

If the solar flux changes or is measured more accurately in the future, 
our method requires simply to modify C 0 to 

Cq 1 = C 0 + log I' - log 1 (7) 

where I and I f are the old and new solar flux, respectively, whereas all 
other coefficients (i =£ 0) remain unchanged. The correction can also be 

made, if necessary, for the effects of changes in the sun-earth distance; 
the earth is about 3.28 percent closer to the sun in early January than in 
early July. The same method can also be used to calculate atmospheric heat- 
ing by diffuse radiation reflected or scattered from clouds or ground. Using 
a known (or assumed) albedo, I T is then calculated from the downward radi- 
ation flux arriving at the reflection heights. The angular integration for 
the diffuse radiation is taken into account by using an air mass factor of 
1.66. Using the more recent value 1.9 by Lacis and Hansen (ref. 8 ) would 
increase the calculated heating rate by diffuse radiation by ~5 percent. 

Our result is compared with the results calculated by different formulae 
developed by other investigators in figure 2. Since the results are so close, 
the illustration of comparison of each model with our result is shown by 
shifting vertically with each other. 

Manabe and Strickler (ref. 9) have used an analytical equation for 
absorptivity; multiplying the total incident solar energy (~2 cal cm -2 min ” 1 
or ~ 1 . 39 xl 0 6 erg cm “2 sec- 1 ) to the absorptivity would give the total 
absorbed energy. As is seen in figure 2 their result is in good agreement 
with the present study for the range of u between 10” 1 and 1 cm, which 
covers the main part of the stratosphere but gives much smaller absorption at 
larger values of u (^2 cm) and at smaller values of u (<10 “ 3 cm) . Thus 
their equation would cause serious errors and should not be used for the 
region above ~50 to 70 km (the actual height depends upon seasons and lati- 
tudes) , where the path of the solar beam is smaller than ~ 10~ 3 cm, and for the 
region below ~25 km in winter at high latitudes, where the beam path is 
greater than ~2 cm. 

Lacis and Hansen (ref. 8 ) used analytical formulae for absorptivities 
representing the spectrum of ultraviolet and visible radiation. As seen in 
figure 2 , the total absorbed energy calculated for the two spectrum ranges 
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by their formulae is slightly smaller than the corresponding value calculated 
by our formula in most of the range of u except for the range between about 
1 and 10 cm. 

Lindzen and Will (ref. 10) have formulated the differential heating rate 
for the three bands into three different formulae. The absorbed energy can 
be calculated by integrating these equations from 0 to u. Using equations 
similar to those of Lindzen and Will, Strobel (ref. 11) has determined the 
new coefficients to calculate the heating rate for the three bands. The 
total absorbed energy calculated by these equations is compared with our 
result (fig. 2). Both Strobel and Lindzen formulae give absorbed energies 
in good agreement with the result of the present study for the range of u 
between 10“ 2 and 10 cm, which covers most regions of the stratosphere. Our 
formula gives, however, generally larger absorption than most other formulae, 
and part of the reason for this difference may be caused by the fact that our 
analyses include wider spectra, although the ranges are not always clearly 
stated in some models. 


Heating by Water Vapor 


Water vapor is the main absorber of solar radiation in the troposphere. 
Howard et al . (ref. 2) have formulated their experimental results for 
integrated absorption function for seven bands (6.3 y, 3.2 y, 2.7 y, 1.87 y, 
1.38 y, 1.1 y, and 0.94 u) into the following formulae 



A dv = cu 

v 


(P + e)' 


( 8 ) 


for weaker absorption, and 

J A dv - C + D log u + K log (p + e) (9) 

Av v 

for stronger absorption, where u is the water vapor column density in 
g cm” 2 , p and e represent the total pressure and the partial pressure of 
water vapor in mm Hg, respectively. The constants c, d, k, C, D and K are 
given for each band (see Burch et al . , ref. 12, for the revised values for 
these constants) . The absorbed solar energies calculated using these equa- 
tions are shown for each band and for the total of seven bands in figure 3. 
Sudden jumps on some curves are caused by the transition between the range 
where the weaker absorption form (8) is applicable and the range where the 
stronger absorption form (9) is applicable. There is no apparent disconti- 
nuity on the curve for the total of seven bands. We have fitted the total 
curve with equation (5) for the range of u between 10 -7 and 10 2 g cm" 2 ; the 
determined coefficients for the seventh-order polynomial are tabulated in 
table 2. The variation calculated by this equation agrees well with the 
original variation which is shown by a dot-dashed curve in figure 3. 


Many attempts have been made in the literature to express the solar 
energy absorbed by water vapor in simple analytical forms, and the results 
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of calculations by these formulae are compared in figure 4. The formula of 
Miigge and Moller (ref. 13) is essentially the first-order polynomial of the 
type of equation (5) , and changing the units from cal cm" 2 min” 1 to 
erg cm" 2 sec” 1 , their equation can be written as 

log S = 5.079 + 0.303 log mu (10a) 

where m is the air mass factor. Our corresponding equation is expressed as 

log S = 5.007 + 0.429 log mu (10b) 

Korb et al. (ref. 14) have improved equation (10a) and obtained the third- 
order polynomial 

log S = 5.1036 + 0.347 log mu - 0.056 (log mu) 2 - 0.006 (log mu) 3 (11a) 

which can be compared with our third-order polynomial 

log S = 5.0791 + 0.3150 log mu - 0.0416 (log mu) 2 - 0.0029 (log mu) 3 M1 , N 

(lib) 


The result from the Korb formula gives the closest variation to the 
present result, which is very close to the one calculated from the formula 
of Howard et al. As seen from table 3, the mean error of our third-order 
polynomial is about 2.4 times larger than the mean error of the seventh- 
order polynomial; therefore, it is preferable to use the seventh-order poly- 
nomial . 

The analytical forms of Sasamori et al. (ref. 15) and Lacis and Hansen 
(ref. 8) are both based on absorptivities evaluated with the radiation chart 
of Yamamoto (ref. 16); thus, the general tendency of the variation in these 
two results is similar over the range of mu = 10~ 2 - 10 cm for which the 
formulae may have been approximated. Extension of the formulae beyond this 
range causes large discrepancies between the two results. It is seen in 
figure 4 that the formula of Howard et al. (ref. 2) gives smaller absorption 
than most other formulae for the range of u greater than ~10” 2 g cm -2 , 
whereas it gives larger absorption for the smaller u than this value. 


Heating by Carbon Dioxide 

Howard et al. (ref. 1) have used the same formulae as equations (8) 
and (9) to express their experimental data on the absorption by CO 2 in 
eight bands (14.7 y, 5.2 y , 4.8 y, 4.3 y, 2.7 y, 2.0 y, 1.6 y, and 1.4 y). 

The calculated absorption in each band and the total of eight bands are 
shown in figure 5. The large jump appearing in the curve of the 4.3 y band 
causes a discontinuity for the curve of total of eight bands, since the 4.3 y 
band gives the largest contribution. We have smoothed in the region of the 
discontinuity and have then applied the curve fitting process using equa- 
tion (5) for the range of u between 10” 7 and 10 3 cm. The coefficients deter- 
mined for the seventh-order polynomial are tabulated in table 2 along with the 
mean errors in each case of third to seventh-order polynomials. 
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DISTRIBUTION OF GASES 


To calculate the atmospheric heating rate using the analytical formula 
developed in previous sections, we need the data of global distributions of 
absorbing gases. 


Ozone 

The general latitudinal and seasonal variations of the vertical ozone 
distributions up to ~35 km are reasonably well known. We have used data 
analyzed and compiled by Wilcox et al. (ref. 17) based on ozonesonde data 
obtained at various stations in the world from 1957 to 1972. Wilcox et al. 
give the ozone densities for every 2.5 km up to the height of 32.5 km, and 
we have interpolated these values to obtain the ozone densities at every km 
interval. For the region above 30 km up to ~55 km, the satellite measurements 
of backscatter ultraviolet (BUV) solar radiation from the atmosphere can give 
the information of ozone profiles on a global basis. Taking the data at 50 km 
obtained by the Nimbus-4 satellite experiments (I. Eber stein, private communi- 
cation, 1977) and combining it with the ozonesonde data at 32.5 km, we have 
calculated by interpolation and extrapolation the ozone densities at every km 
up to the height of 80 km, assuming the exponential decrease with height. 

The result is shown in figure 6(a) for March (northern spring and southern 
fall) and in figure 6(b) for June (northern summer and southern winter). 


Water Vapor 

The distribution of water vapor mixing ratio in the troposphere is 
expressed in the analytical form of 


p h 2 0 ■ f(z)e ' SX <12) 

where x is the latitude, z is the height and s is the parameter changing 
sinusoidally with season from 0.016 in summer to 0.032 in winter. These 
values of s and the functional form of f(z) are determined based on data 
published by Oort and Rasmusson (ref. 18). The adopted values of f(z) are 
given in table 3. In the stratosphere the mixing ratio is assumed to be con- 
stant at 3.5x10 6 . Since the tropopause is low (~8 km) at high latitudes, 
the above procedure would cause sudden large changes in the water vapor mix- 
ing ratio at ~8 km at latitudes above ~60°, as is seen in figures 7(a) and (b) . 


Carbon Dioxide 

We assume the constant mixing ratio of 3.2xl0“ 4 for C02* Recent model 
calculations (Shimazaki, unpublished result) indicate that the mixing ratio 
decreases slightly with height (from 3.2xl0~ 4 at the ground to ~1. 86X10 - ** 4 at 
60 km) , but this slight change should not greatly affect the calculation of 
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the atmospheric heating rate, since the C0 2 absorption is not the primary 
source of the atmospheric heating. 


RESULT AND DISCUSSION 


Calculating the absorbed energy in the atmosphere by equation (5), we 
have taken into account the dependence of absorption on pressure by introduc- 
ing the effective mass of the radiative absorbent 


u 


eff 



(13) 


where P Q is the normal pressure of one atmosphere (1013.25 mb). For H 2 0 
and C0 2 the power n can be evaluated from k/d in equation ( 8 ) and from 
K/D in equation (9) , and the average values for all bands under consideration 
give ~0.6 for H 2 0 and ~0.8 for C0 2 . The power n for 0 3 is taken as 0.2 
referring to Kondratyev (ref* 19). 

The global average of height profiles of the heating rate calculated for 
each constituent for March is illustrated in figure 8 . The figure demon- 
strates the strong effect of n on the heating rate by 0 3 absorption in 
regions above ~40 km. The value of 0.4 has been obtained for n by 
McClatchy et al. (ref. 20), but it is inappropriate to use this value for the 
ultraviolet range, since it would give too small heating rate above ~35 km as 
is seen in figure 8 . McClatchy 1 s value is based on the experimental data 
on ozone transmittance for the infrared radiatidn. The result for n = 0 
gives too much heating above ^50 km particularly at high latitudes. 

Heating by the reflected solar radiation is calculated assuming the 
global average of the effective albedo to be 0.25. The effective albedo 
takes into account the combined effect of reflection from clouds and the 
ground surface and is useful in calculating the heating by the reflected 
waves above the clouds. The present study is concerned mainly with the 
stratosphere and mesosphere, and the concept of the effective albedo should 
give the reasonable heating rate in these regions. It is essential, however, 
for the troposphere to take into account the cloud albedo and the surface 
albedo separately and to calculate the radiative heat fluxes between the cloud 
and the ground surface and also within the clouds. 

The calculated heating rates by reflected waves due to H 2 0 and C0 2 
absorption are very small in the stratosphere (fig. 8 ). Furthermore, since 
the albedo of clouds is relatively small in the infrared and also since 
water vapor droplets in clouds absorb radiation almost at the same spectrum 
as water vapor, the heating by reflected waves due to H 2 0 and C0 2 absorption 
should be smaller than those shown in figure 8 and can be entirely neglected. 
On the other hand, the albedo of clouds is high in the visible and ultra- 
violet, therefore heating by the reflected waves due to the O 3 absorption 
cannot be neglected. In fact, as is seen in figure 8 , this heating is com- 
parable to the heating by C0 2 (direct radiation) in the lower stratosphere, 
and is about 22 percent of the O 3 heating by direct radiation at 25 km. 
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Applicability of our formula is limited to the range of u greater than 
10 “ 7 g cm “ 2 for H 2 O and 10 ” 7 cm atm NTP for O 3 and CO 2 ; these critical values 
are met somewhere above ~40 km for H 2 O and above ^80 km for CO 2 and O 3 (see 
the top abscissa of figs. 1, 3, and 5), although these heights depend upon 
the season and latitude and should be much higher near the sunrise and sunset 
times. The heating rate in the regions outside of these critical heights are 
not shown in figure 8 and are neglected in our calculation of atmospheric 
heating. Fortunately, however, contributions from these heatings are gener- 
ally small and their neglect should not affect the global distribution of the 
total heating rate in the region below 80 km. The upper limits of the for- 
mula’s applicability are 10 2 g cm ~ 2 for H 2 O, 10 2 cm atm NTP for O 3 and 10 3 cm 
atm NTP for CO 2 . The column densities of the absorbents are usually within 
these limits. 

The heating rate due to CO 2 and H 2 O absorption represented in terms of 
°K per day is almost constant above certain heights (fig. 8 ). This phenom- 
enon is caused by the assumption of the constant mixing ratio, that is, these 
constituents’ concentrations change with height in proportion to the mass 
density of atmosphere. Since the heat capacity also decreases with height in 
proportion to the mass density of the atmosphere, the rate of change in tem- 
perature should remain almost constant with height, although the actual 
absorbed energy decreases with height. 

Latitudinal distributions of the diurnal heating rate are illustrated in 
figures 9(a) and 9(b) for March and June, respectively, in terms of tempera- 
ture change per day. The top 1 to 2 km region is affected by the assumption 
on the total amount of 0 3 above the upper boundary: we assume that the O 3 

density decreases above 80 km with the scale height of its own mass. The 
actual distribution is complicated and has some bulge around 80 km (see, for 
instance, Schimazaki and Laird, ref. 21). For March, the global maximum 
occurs around 40 to 45 km in the equatorial region. Whereas the distribution 
is almost symmetrical around the equator in regions below this height, a 
marked asymmetry appears in higher regions; the larger heating rate is cal- 
culated at high latitudes in the southern hemisphere. This is certainly 
related to the larger ozone density calculated in that same area (see 
fig. 6 (a)). 

It is well known that there is a worldwide maximum of the total ozone 
density at the northern high latitudes in March and April. This tendency is 
clearly seen in the ozone density at ~30 km in figure 6 (a); that is, the 
ozone density observed at the northern high latitudes is appreciably larger 
than the corresponding values at the southern high latitudes. It is not cer- 
tain, however, whether or not this tendency should remain the same in the 
higher regions since contribution of that region to the total ozone density 
is very small. 

The ozone density above 50 km shown in figure 6 (a) is calculated by 
extrapolation from the ozonesonde value at 32.5 km and the BUV value at 
50 km. The larger ozone density above 50 km at the southern high latitudes 
results from the fact that the ozone density is relatively low at 32.5 km, 
whereas the value at 50 km is about the same between the southern and 
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northern high latitudes. There may be little hemispheric difference for 
ozone density in the entire region above 50 km, but the larger ozone density 
above 50 km calculated in the southern high latitudes is consistent with the 
smaller density in the lower region, since the larger ozone density at higher 
region should absorb more radiation, making the radiation available for ozone 
production smaller in the lower region. However, it is true that there is a 
great uncertainty on the latitudinal distribution of the ozone density (and 
therefore the heating rate) above 50 km. 

Calculating the heating rate assuming the uniform distribution of ozone 
density with latitudes represented, for instance, by the distribution at 
35° N, we obtain the uniform distribution of the heating rate with latitude 
above ~50 km. The distribution below this height is affected very little; 
the general features are similar to those shown in figure 9(a). It is evi- 
dent that the global distribution of heating rate is sensitive to the ozone 
density distribution. Unfortunately, there is no direct measurement of global 
distribution above 50 km. Rocket observations provide scant data for limited 
locations and time period. Since the atmospheric scattering is very weak 
above ~55 km, there is a practical upper limit of the altitude from, which 
ozone information may be obtained by the BUV method. A good theoretical 
model would be useful in filling this information gap. 

Heating below ~20 km is mainly due to H 2 O absorption, and the high 
density of H 2 0 near the surface in the lower latitudes shown in figure 7(a) 
causes an increase of the heating rate in that region. 

For June there is a global maximum of the heating rate at ~50 km near 
the north (summer) pole (see fig. 9(b)). Above this height, we again see 
the effect of the nonuniform distribution of the ozone density shown in 
figure 7(b). Below ~20 km where the heating by H 2 0 is dominant, we see 
a large heating in the summer at high latitudes, particularly near the 
tropopause (~8 km) . This is caused by the sudden change of the humidity 
from the troposphere to the stratosphere. Since the tropopause is relatively 
low at high latitudes, and since we have assumed the same small mixing ratio 
for H 2 0 for the entire stratosphere, the model has produced a very large 
sudden change of the H 2 0 concentration at the tropopause at high latitudes 
(see fig. 7(a) and (b)). 

The horizontal distribution of the heating rate is shown in figure 10(a) 
through (h) at every 10 km for the month of June. At the surface, the maxi- 
mum heating occurs at ^10° north of the equator (summer) around noon (180° 
of longitude) and decreases rather steadily in all directions. At 10 km, a 
slow variation occurs from the equator to ~55° N during the midday period, 
and a sudden decrease appears at the northern end. This is caused by the 
sudden change of H 2 0 concentration; the 10 km level is in the troposphere 
(wet atmosphere) at the lower latitudes, whereas it is in the stratosphere 
(dry atmosphere) at the higher latitudes. 

Above 20 km, the ozone heating is the dominant heat source, and the 
single maximum appears in the distribution at 30 and 40 km. There are two 
maxima, however, one larger maximum at the northern summer high latitudes 
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("'■'50 to 70°) and a smaller maximum at the southern winter latitudes 
C — 35 to 50°) at 20 and 50 to 70 km. The latitude where the maximum appears 
shifts towards higher latitudes as the altitude goes up. These features of 
the maximum heating rates directly reflect the similar features of the single 
and double maxima in the distribution of the ozone density (see fig. 6(b)). 


CONCLUDING REMARKS 


This study presents a convenient analytical formula for calculations of 
the atmospheric heating rate. Its form is identical to that used for cal- 
culations of the photodissociation rate at the 0 2 Schumann- Rung e band 
system (Shimazaki et al., ref. 3); it should be useful for reducing the com- 
putational time and memory allocation for computer programs in heating and 
dissociation calculations in three-dimensional global circulation models. 

We plan to work out similar formulation for dissociation calculations for 
many other molecules at different wavelength ranges. 

The mean relative error of our formula for calculations of ozone heating 
is ~2.4 percent of the exact calculations, and our formulae for H 2 0 and C0 2 
are within 1% of the calculations by formulae of Howard et al. (refs. 1 and 
2). The formulae of Howard et al. can predict their observations to ±3%. 

The larger error for ozone heating calculations than for H 2 0 and C0 2 heating 
calculations is caused mainly by the fact that we calculate the heating rate 
due to the three bands by a single equation. These three bands are at very 
different wavelength ranges. Tables 1(a) and 1(b) include the coefficients 
that should be used for each band separately. Using such coefficients can 
reduce the mean error for ozone heating calculations; it would also be useful 
for some specific purpose, such as detailed calculations of the change in the 
heating rate at the ozone cut-off range near 3000 A, when the ozone density 
changes largely from the currently observed value. 

Our formula can be used only for the range of column density greater 
than 10” 7 cm atm NTP for 0 3 and C0 2 , and 10” 7 g cm -2 for H 2 0. These values 
cover most regions below 80 km except for the H 2 0 heating above "~40 km. 
However, the heating by H 2 0 in this region should not affect the total 
heating rate. The ozone heating by the reflected solar short waves contrib- 
utes about 20 to 30% of the heating by the direct solar radiation in regions 
of 20 to 30 km. 

The calculated atmospheric heating rate and its global distribution are 
reasonable, but they are sensitive to the ozone density distribution in the 
stratosphere and the mesosphere and to the water vapor distribution below 
20 km. The global ozone distribution above 50 km is not known well, and 
observational as well as theoretical studies for this region should be 
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important. The distribution of water vapor in the troposphere and its change 
through the tropopause are also important for calculating the realistic 
distribution of the heating rate below ~20 km. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif. 94035, July 19, 1978 
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TABLE 1(a).- COEFFICIENTS OF THE SEVENTH-ORDER POLYNOMIAL FOR HEATING CALCULA- 
TIONS DUE TO 0 3 ABSORPTION AND MEAN RELATIVE ERROR FOR THIRD- TO 
SEVENTH-ORDER POLYNOMIALS (CGS UNITS) 


Spectrum 

range 

Hartley band, 
2400-3000 A 

Huggins band, 
3000-3600 A 

Chappuis band, 
4000-8500 A 

Total, 
2400-8500 A 


—9 — 1 \ 

Incident solar energy (erg • cm z sec ) 


I 

1.406 (4) a 

5.915(4) 

6.918(5) 

8.097(5) 



Coefficients 



Co 

4.14292071(0) 

4.23361218(0) 

4.43346790(0) 

4.75812947(0) 

C 1 

-5 . 33910483 (-2) 

4.12314108 (-1) 

9. 50434910(-1) 

4. 93805176 (-1) 

c 2 

1. 38771166 (-2) 

-1. 60939373 (-1) 

-6 . 23836459 (-2) 

1 . 26465765 (-1) 

C 3 

5.36079604 (-2) 

1. 26356978 (-2) 

-3.20014961(-2) 

2 . 10425653 (-2) 

c 4 

-8. 90746861 (-3) 

9. 98975983(-3) 

-7.37424271(-3) 

-2.45982304(-2) 

C 5 

-6.07296255 (-3) 

6 . 52086922 (-4) 

-6 . 81798995 (-4) 

-7 . 96267282 (-3) 

C 6 

-8.42668609 (-4) 

-1.25888537 (-4) 

-1.89154394(-6) 

-8. 71717239 (-4) 

c 7 

-3. 75005646 (-5) 

-1.27022887 (-5) 

2.17012054(-6) 

-3.24914714(-5) 

Mean relative errors 

E(3) 

1. 264 (-1) 

4.413 (-2) 

5. 796 (-2) 

1.043(-1) 

E(4) 

5. 294 (-2) 

4.429(-2) 

2. 708 (-2) 

6 . 865 (-2) 

E(5) 

4.339 (-2) 

2. 348 (-2) 

9. 653(-3) 

6.588(-2) 

E(6) 

2. 365 (-2) 

1.361 (-2) 

5 .576 (-3) 

2 . 795 (-2) 

E(7) 

1. 678(-2) 

1.255 (-2) 

5.583(-3) 

2.419 (-2) 


ct 

A(b) should read as A x 10 


b 



r 


TABLE 1(b).- COEFFICIENT OF THE SEVENTH-ORDER POLYNOMIAL FOR HEATING 
CALCULATIONS DUE TO 0 3 ABSORPTION (MRS UNITS) 


Spectrum 

Hartley band. 

Huggins band. 

Chappuis band. 

Total, 

range 

2400-3000 A 

3000-3600 A 

4000-8500 A 

2400-8500 A 


Incident solar energy (J 

♦ • sec * ) 


I 

1.406(1)“ 

5.915(1) 

6.918(2) 

8.097(2) 



Coefficients 



C 0 

1.12492560(0) 

1.68658866(0) 

2.68914243(0) 

2.71161737(0) 

Cl 

-3. 04055682 (-1) 

2.62163034(-1) 

2 . 69717248 (-2) 

-3.53904066 (-1) 

c 2 

-5.91532217 (-1) 

1.68046817 (-1) 

-4.84914017 (-1) 

-1.20569660(0) 

C 3 

-4.16397584 (-1) 

9 . 13818058 (-2) 

-1 . 17354777 (-1) 

-6 . 51920173 (-1) 

c 4 

-1. 30697369 (-1) 

5. 40067599(-3) 

-1. 36980916(-2) 

-1. 65625605 (-1) 

C 5 

-1.93350333 (-2) 

-1.92556778 (-3) 

-5. 22207398(-4) 

-2.11525633(-2) 

C 6 

-1. 36767651 (-3) 

-3. 03720579 (-4) 

2. 84901436(-5) 

-1. 3265 9784 (-3) 

. C 7 

-3. 75005646(-5) 

-1. 27022887 (-5) 

2.17012054(-6) 

-3. 24914714 (-5) 


a A(b) should read as A x 10^. 
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TABLE 2.- COEFFICIENTS OF THE SEVENTH-ORDER POLYNOMIAL FOR HEATING CALCULA- 
TIONS DUE TO H 2 0 AND CO 2 ABSORPTION AND MEAN RELATIVE ERROR FOR THE 
THIRD- TO SEVENTH-ORDER POLYNOMIALS 



CGS Units 

MKS Units 

Spectrum 

h 2 o 

C0 2 

H 2 0 

co 2 

range 







Incident solar incident 


I 

3.292(5)“ 

1.162(5) 

3.292(2) 

1.162(2) 



Coefficients 



Co 

5.05447794(0) 

3.53768885(0) 

1.73903147(0) 

1.10124356(0) 

Cl 

2.91792812 (-1) 

3. 70827458(-l) 

3.53646788 (-1) 

2 . 46467496 (-1) 

c 2 

-8. 46611844 (-3) 

-5. 91767097 (-2) 

-5 . 04237590 (-2) 

1. 95386476 (-2) 

c 3 

1. 37174940 (-2) 

9. 35607274 (-4) 

8. 61937376 (-3) 

2. 20605560(-2) 

c 4 

-4. 03685379 (-3) 

3.04430731 (-3) 

4.5107195K-3) 

4 . 65717546 (-4) 

C 5 

-3.04986573(-3) 

2. 00068263 (-4) 

-5. 59632427 (-4) 

-8. 62659134(-4) 

C 6 

-5. 10273302 (-4) 

-5.18423869(-5) 

-3.19804466 (-4) 

-1. 25278846 (-4) 

C 7 

-2. 72098338(-5) 

-5. 24546137 (-6) 

-2 . 72098338 (-5) 

-5. 24546137 (-6) 



Mean relative errors 


E(3) 

1.827 (-2) 

2. 780(-2) 

1.827 (-2) 

2. 780(-2) 

E(4) 

1. 771 (-2) 

2.613 (-2) 

1 . 771 (-2) 

2. 613(-2) 

E(5) 

1.560 (-2) 

1.238 (-2) 

1.560(-2) 

1. 238 (-2) 

E(6) 

1.443 (-2) 

9.552(-3) 

1 . 443 (-2) 

9 . 552 (-3) 

E(7) 

7.681 (-3) 

8. 329 (-3) 

7 . 681 (-3) 

8.329 (-3) 


a A(b) should read as A x 10^. 
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TABLE 3.- ASSUMED H 2 0 MIXING RATIO AT EQUATOR 
(FUNCTION f(z) IN EQ. (12)) 


Z (km) 

f 

(z 

) 

0 

1.09 

X 

10~ 2 

1 

6.84 

X 

10 -3 

2 

4.65 

X 

10~ 3 

3 

3.42 

X 

10 -3 

4 

2.33 

X 

10 -3 

5 

1.39 

X 

10 -3 

6 

8.55 

X 

10 -4 

7 

5.44 

X 

10 -4 

8 

3.11 

X 

10 -4 

9 

1.55 

X 

10 -4 

10 

9.00 

X 

10~ 5 

11 

5.3 

X 

10~ 5 

12 

3.05 

X 

10 -5 

13 

1.75 

X 

10 -5 

14 

1.0 

X 

10 -5 

15 

5.8 

X 

10~ 6 

16 

3.5 

X 

10~6 
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Figure 1.- Absorption by o zone at Hartley band, Huggins band, and Chappuis 
band. The total of the three bands is compared with the result calculated 
by our formula. The abscissa indicates the ozone column density, and the 
corresponding height scales are shown at the top for three representative 
seasons and latitudes. 
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OZONE COLUMN DENSITY, cm-atm NTP 


Figure 2.- Comparison of ozone absorption calculated by our formula with those 
calculated by various other formulae in the literature. Arrows indicate 
the ordinate at 10 2 erg* cm” 2 sec” 1 for each pair of curves. 



HEIGHT, km 


40 

30 

20 15 10 & 

5 

0 WINTER, 65° 

1 

30 

l 

20 

I l I 

16 '15 10 

I 

5 

1 

0 EQUINOXES, 0 C 


I f I l I I 1 


40 30 20 15 10 & 5 0 SUMMER, 65° 

i i i \ i I 1 



Figure 3.- Absorption by water vapor at seven bands in the range of 0.94 to 
6.3 y calculated by the formulae of Howard et al. (ref. 2). The total of 
seven bands is compared with the result calculated by our formula. 
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H : HOWARD-BIRCH-WILLIAMS (1956 b) 

M : MUGGE AND MOLLER (1932) 

K: KORB et al. (1956) 

S : SASAMORI etal.(1972). 

L : LACIS AND HANSEN (1974) 



Figure 4.- Comparison of water vapor absorption .calculated by various formulae 
in the literature. Our formula gives the result very close to the curve 
of Howard et al. 
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Figure 5.- Absorption by carbon dioxide at eight bands in the range of 1.4 to 
14.7 y calculated by the formulae of Howard et al. (ref. 1). The total of 
eight bands is compared with the result calculated by our formula. The 
area of discontinuity at u ^ lO” 1 cm is smoothed out before the curve- 
fitting process is applied. 
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(a) March. 


h 2 o number density 



(b) June. 

Figure 7.— Latitudinal distribution of the water vapor density 



GLOBAL AVERAGE, MARCH 



HEATING RATE, K day -1 


Figure 8 .- Height variation of the heating rate due to O 3 , H 2 0 and C0 2 
absorption, representing the global average for March- The heating 
rate due to reflected waves is calculated assuming an effective albedo 
of 0.25. A very large pressure effect appears on the ozone heating 
above ~40 km. The value of n = 0.2 was adopted in the present study 









HEATING RATE, erg/cm 3 • sec NORTHERN SUMMER 



I 30 60 90 120 150 180 210 240 270 300 330 360 

LONGITUDE, deg 


(c) 20 km 


HEATING RATE, erg/cm 3 • sec NORTHERN SUMMER 



(d) 30 km 


Figure 10.- Continued. 
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Figure 10.- Concluded. 
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